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Abstract 

We consider a holographic model of chiral superfluidity whose bulk is Einstein Yang- 
Mills and compute viscosity and conductivity responses away from the probe limit. We 
calculate Hall viscosity and analyze its relationship to the superfluid density and the 
shift. We find that the relationship between these quantities derived from effective 
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equal to ^1/2. At low temperatures the system develops a Lifshitz throat, indicating 
an anisotropic scaling symmetry in the infrared dynamics. 
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1 Introduction 


The promise of providing a framework to construct Majoranna fermions and its potential ap¬ 
plication to fault tolerant quantum computing lias made high temperature chiral superfluids 
an attractive subject of study in recent years [1, 2, 3]. While traditional superconductiv¬ 
ity is well understood via the BCS pairing mechanism, high temperature superconductivity 
cannot be explained by conventional electron-phonon interactions that dominate the pairing 
mechanism in BCS. 

There has been much advancement in the understanding of these systems [4, 5, 6], how¬ 
ever a universal understanding which can explain different physical aspects is not currently 
in hand. As an example, the stress response at finite temperature has not been fully ex¬ 
plored. Given the possible practical application of these systems in quantum computing, it 
is important that their tensile properties be fully analyzed and understood. 

One possible venue for arriving at universal predictions is by considering the effective field 
theory (EFT) that governs the dynamics at very low temperature and energies. The main 
guiding principle of this approach is the underlying symmetries of the system. For example, if 
we consider consider a superfluid system controlled mainly by a quantum critical point with 
Lorentz invariance, the EFT must also be Lorentz invariant. These symmetries are often not 
powerful enough to completely constrain the low energy dynamics of the system, however, 
in some cases when combined with topological data they can provide predictions that are 
either universal or can relate various transport coefficients. This approach was applied to 
superfluid systems at zero temperature in [7, 8, 9]. In particular it was shown that the ratio 
of hall viscosity to the superfluid density can be determined by topological considerations, 
i.e. it is given by the superfluid shift, which relates the ground state vortex number of the 
system to the spatial topology. Specifically, the ratio of Hall viscosity to superfluid density 
is equal to =f 1/2 f° r p±ip superfluid. However, by their very essence, these analyses are not 
expected to be valid at finite temperature and since they rely on an action principle, they 
also cannot describe any dissipative physics. 

Traditionally, in areas of condensed matter physics where a careful first principle anal¬ 
ysis is not available or difficult, holography has stepped in to provide valuable intuition by 
providing access to a very different set of universality classes than that of weakly interact¬ 
ing quasiparticles. Since these systems are the result of deforming conformal theories, the 
inherited conformal symmetry further constrains the transport properties of the system. 

In this paper we use a holographic model to compute the stress response of relativistic 
chiral superfluids. The model we consider is Einstein-Yang-Mills with group SU( 2) [10, 11], 
where, after breaking SU( 2) —>■ U( 1) the W-bosons in the bulk provide a charged vector 
operator in the boundary. Studying stress response requires looking at the system away 
from the probe limit (where the charged and neutral sectors decouple) and so we must allow 
for backreaction. We also analyze the electronic responses that were previously computed 
in the probe limit and discuss how they change. The main result of this analysis is that the 
ratio of Hall viscosity to superfluid density derived in the EFTs persists at finite temperature 
all the way to Tq- 

The outline of the paper is as follows. In section 2, we review the main results of the 
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effective field theory that governs the dynamics of the system at zero temperature. In section 
3, we discuss the background geometry of the system and analyze the behavior of various 
background quantities. At low temperatures the system develops a new scaling symmetry 
with nontrivial dynamical critical exponent z. In section 4, we calculate Hall and shear 
viscosities and conductivities. We conclude in section 5. 


2 Review of effective field theory predictions 


In this section we review some of the results from the effective theory description of chiral 
superfluids that holds at zero temperature [8]. We work in 2 + 1 dimensions. Defining the 
dual gauge held a M as the hodge dual of the superfluid current j — A */, the action up to 
first order in derivative expansion 1 is: 


S 



-e{n) + 


1 
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+ C (nje^PupdvUp + nap , 


( 2 . 1 ) 


where Ap is the electromagnetic potential, the number density n is defined as n 2 = —j^jp = 
+/ 2 /87t 2 , e{n) gives the energy density and £(n) is an undetermined function of the number 
density. Finally, the current is the Euler current[12, 8]: 


1 

87T 


= _ £ ^P £ *Pl u 


7 pUry 


:R 


■isp(3 , y 


( 2 . 2 ) 


where is the normalized direction of the superfluid current = nu A We note that 
the total charge associated to the Euler current on a spatial section is a half of the Euler 
characteristic of the surface [8]. 

It is important to note that since a M is the dual of the superfluid current, its charge 
counts the vortex number in the superfluid. By taking a variation of the action (2.1) with 
respect to a o, we can see that there is a relationship between the vortex number, the total 
flux and the Euler characteristic of the space: 

N V = N' P + /c|. (2.3) 

This is in fact a constraint built into the effective action from topological considerations [13, 8] 
and in the case of a p ± ip superfluid, it can be summarized as the necessary relaionship 
between the flux and the topology such that the vector order parameter can be smoothly 
defined on the entire manifold. In terms of the effective theory parameters, this means: 


p ± ip : k = =f2. 


(2.4) 


Defining the superfluid shift S as the net flux required to have no vortices when the superfluid 
is placed on a sphere, we see that 2 S = =f 2 for p ± ip. 

1 In this power counting scheme is taken as 0(p~ 1 ) because n is O(p 0 ). The background fields Ap and 
gp V are also taken as O(p 0 ). 

2 Comparing to [8], we have put the superfluid charge equal to 1. 
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Of particular interest are the predictions of this effective action with regards to the Hall 
conductivity and viscosity at zero frequency. For the p + ip case, 

n 

TlH = - r 

nq 2 2 q 2 ((n) 

axy = ~~i^ + 

While the prediction for the Hall viscosity is universal, the Hall conductivity depends on 
the non-universal function ((n). Requiring the action to be Weyl invariant fixes the this 
function to be linear ((x) = (qx and we get: 

a xy = { 2Co-l)^. (2.7) 

/i 

In what follows, we verify these predictions in the holographic setup and find that the Hall 
viscosity result is in fact robust at finite temperature. 


(2.5) 

( 2 . 6 ) 


3 Holographic chiral superfluids 

In this section we review the standard construction of holographic chiral superfluids [10, 11]. 
The bulk theory is Einstein-Yang-Mills with group SU( 2), 3 where by turning on a chemical 
potential for a generator r 3 we explicitly break the global symmetry group of the boundary 
theory to U( 1). The W-bosons in the bulk provide us with natural charged vector operators 
in the boundary theory, which we use to construct a p + ip order parameter spontaneously 
breaking the remaining U( 1). This spontaneous symmetry breaking is vector-like and breaks 
parity and time reversal. 


3.1 The superconducting phase 

Consider Einstein-Yang-Mills theory in asymptotically AdS^ space: 


S = J d 4 x^/^g 


1 

2k 2 




l 

W 


(3.1) 


where + f^A^Al, and we chose G = SU( 2), f abc = e ( 

Einstein’s equations are simply 


abc 


2k 2 


— ( R,, - -Rg, v - —g iW ) = T M „, T, v = — [ - -F a pa F a ^ ) , (3.2) 




ri a X 
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and the Yang-Mills equation is 


^ ^paiiis _j_ j?abcj^b ecu _ q 


(3.3) 


3 Note that even though this holographic p + ip superfluid is unstable to a rotation breaking p x phase 

[14, 11], this is not an problem when calculating the response functions under consideration here. 
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In what follows it will be advantageous to work with complex coordinates and generators, 


x + iy , 

2 =- r ^ = 


T 1 ± IT 2 


y/2 ’ y/2 ' 

Our normal phase is the standard Reissner-Nordstrom black hole, 

r 2 rH{2g 2 r 2 H + k 2 /i 2 L 2 ) r^K 2 ^ 2 


(3.4) 


ds z = -fdt 2 + -j-+2r 2 dzdz, f = 


L 2 


2 g 2 L 2 r 


+ " 77 , A = /i(l -r H /r)r 3 dt. 

(3.5) 

In this normal phase, the non-normalizable chemical potential explicitly breaks SU{2) down 
to 17(1). At low enough temperature and large enough charge there may be a dynamical 
instability to developing vector hair. Consider the following ansatz: 


ds 2 = — h{r)dt 2 + 2r 2 dzdz + 


dr 2 

f( r Y 


Px + ipy : A = cj)(r)T 3 dt + w(r)(r dz + r + dz), 
p x - ip y : A = (p(r)r 3 dt + w(r)(r + dz + r~dz). 

This is invariant under a combined rotation and gauge transformation 

p x ±ip v :z^ e"z, U = 


(3.6) 


(3.7) 


The equations of motion for the background fields are independent of chirality and reduce 
to 
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g 2 \2r 3 / 2 h 


w’ = 0 . 


(3.9) 

(3.10) 

(3.11) 


Note that our ansatz has two scaling symmetries [10], 


r —> ar , (£, x, y ) —» (t, x, y)/a, f —> a 2 /, h ci 2 h, w —$■ aw, (j) —» acj), 


t —^ bt , h —>■ 5 2 /i, <f) —>■ bcj). (3.12) 

By considering perturbations of w on the Reissner-Nordstrom background and looking for 
normalizable solutions we can extract the critical temperature as a function of g = gL/n, 
which we show in figure 1. We find that there is a critical value of the coupling gc ~ 1.045 
where the critical temperature vanishes. 
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Figure 1: Plot of the critical temperature as a function of g. 


When considering the system at finite temperature we look for solutions with a regular 
horizon corresponding to / vanishing linearly at r = t'h and we will pick a gauge where 
4>(r H ) = 0 as well. Asymptotically the fields behave as 


„ r 2 M r 2 M 

f — To -+ ■■■, h = — -+ 


f~h = 0(r~ 2 ) 


L 2 r ’ L 2 r 

w = w 0 + Wi/r + ..., (/> — p — p/r + ... 


(3.13) 


Following the standard holographic dictionary we interpret p as the chemical potential and 
Wo corresponds to an external source for the nonabelian current whose expectation value is 

( Jbdy) = + dzr + ). (3.14) 

To have spontaneous U( 1) breaking we consider only solutions with Wq = 0. The entropy, 
temperature, energy density and total charge density are 


5 = 


27rr| 

9 ! 


T = 


Vf'( r n)h'(r H ) 


47r 


M _ p 

£ o" i Ptotal o • 

K z g z 


(3.15) 


We can define the normal fluid charge density via the flux on the horizon and similarly the 
superfluid charge density as the difference between the total and the normal fluid charge 

[14], 

r 2 

Pn V0 (), Ps Ptotal Pn • (3.16) 

r 

Plots of the superfluid density are given in figure 2. As expected, the ratio of superfluid to 
total charge goes to one at zero temperature and vanishes linearly at Tq as in mean held 
theory. 

We should note that defining the superfluid density through the residue of the zero 
frequency pole of various correlators can be misleading. For example the a xx pole only 
agrees with our definition at zero temperature, and even above Tq the pole is nonzero as 
our system is exactly translationally invariant so there is infinite DC conductivity. While 
we strongly suspect that the definition we use for superfluid density will agree with the 
normalization of the Goldstone two-point function, we leave confirming this to future work. 
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Figure 2: Plot of superfluid fraction and superfluid density (in units of vs temperature. 


3.2 Low temperature and the Lifshitz throat 


At very low temperatures the geometry develops a Lifshitz throat. We can most easily see 
this after noticing that at low temperatures the solutions behaves very similarly to an s-wave 
holographic superconductor, where most of the electric charge in the solution is carried in 
the hair outside of the black hole and not on the horizon itself [15]. In the zero temperature 
limit, the degenerate horizon carries no charge and has vanishing area. In the extreme near¬ 
horizon limit, we expect a scaling solution with 0 = 0. Examining the equations of motion 
with 0 = 0 and again using the dimensionless coupling g = gL/n we find the solution 


6 — g 2 (z — l) 2 r 2 


/ = 


2(2 + z) 

where z is the real root of the cubic equation 


L2 , h = h 0 r 2z , w = gy/z - 1-, 


z 3 + z 2 + 9 


A- 3 ( 2+ />= 0 . 
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(3.17) 


(3.18) 


and should not be confused with the complex coordinate on the plane. Note that a positive 
root exists for any real g and z > 1. One can also check that the coefficient in front of / in 
(3.17) is always positive. We plot z(g) in figure 3. This throat has an anisotropic Lifshitz 
scaling symmetry [16], 

t —> A z t, x —> A x,r —> rf A. (3.19) 


From this Lifshitz throat we can turn on a perturbation in 0, which (at zero frequency 
and momentum) can have scaling behavior 

0 ~ ar- (1 -" /2+i/) + br~ < ' 1 ~ z/2 ~ u \ v = y/z (4 + z )/4 + 3. (3.20) 

The a mode diverges at small r and vanishes at large r, while the b mode vanishes at small r 
and diverges at large r. We can therefore turn on a 6-type perturbation and its backreaction 
will flow us to AdS^ at large r. 
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Figure 3: Plot of the dynamical critical exponent z and IR dimension A as a function of g. 


We can also consider relevant perturbations in the other background field modes. There 
are normalizable static perturbations of /, x, w which have the following behavior at large 
radius: 

M h 0 L 2 (2 + z 2 )(3 + 2z + z 2 ) M La/3 + 2z(2 + z) M 

= 6z(2 + z) W = 2VGzVTTir^ rz ' 

Such an asymptotically Lifshitz solution has a mass density 



MLy/h Q y/Z + 2 + z 2 
'/Qk? 


(3.22) 


If we consider an asymptotically Lifshitz black hole at finite temperature, it must have 
/ ~ h ~ (r — r H ), and we can use the scaling symmetry (3.12) to argue that 


S ~ Vh, T ~ r z H , e ~ r£ z , 
which can be repackaged into a more familiar form 

S ~ T 2/z , e ~ t 1+2/z , 


(3.23) 

(3.24) 


as expected for a fixed point with dynamical critical exponent z. The numerical evaluation 
of these quantities is given in plot 4 and matches this scaling. We note that due to numerical 
difficulty, the smallest value of g that we considered was g — 5. 


3.3 Superfluid on a sphere 

The superfluid shift, defined in [8] for relativistic superfluids in 2+1 dimensions, is the net 
flux required to have no vortices when the superfluid is placed on a sphere. For a chiral 
p-wave superfluid with charge equal to 1, the shift is related to the chirality: 


p x ±ip y : S = + 2 . 


( 3 . 25 ) 















Figure 4: Plots of entropy and energy density versus temperature respectively in units of 

T ^ 

and energy density 


As the holographic model under consideration is meant to describe a p ± ip superfluid, 
we may ask whether or not it can reproduce this result. The answer it that it indeed 
does. First, as the arguments for this constraint are purely topological they would straight¬ 
forwardly apply to the holographic setup as well. However, putting this system on a sphere, 
one finds the constraint is also embedded in the bulk equations of motion. To see this, 
we look for vortex-free solutions in asymptotically global AdS. Parts of the calculation are 
simpler in projective coordinates, where the sphere metric is given by: 

Or/ 7r j y 

dQ 2 = d6 2 + sin 2 Odcf) 2 = —^ — v^e^tan [9/2 ]. (3.26) 

(1 + zz/2) 


We also need to turn on a constant magnetic flux through the sphere given (in the upper 
hemisphere) by: 


Ah, 


A D AA(1 — cos 9)dc/)T 3 = « 


. N(j, zdz — zdz 3 


9 , - • - (3-27) 

2 + zz 

We generalize the ansatz (3.6) for the sphere such that it retains the same rotation properties 
about the north pole as flat space: 


p±ip : Ad w(r) (r T dz + r^dz') H(zz). 


(3.28) 


The full ansatz for vortex free solutions is 

ds 2 = — h{r)dt 2 + r 2 dVl 2 + -jy-r-, 

f\T) 

A = t 3 fcj)(r)dt + i — zdz\ + _|_ T ^dz) H(zz)w(r). (3.29) 

V 2j 2j | zz J 

Plugging this into the equations of motion, we find that the bulk Yang-Mills equations 
constrain both and H. First, looking at the (.. r component, we find the constraint 


H\zz) 


N^Hjzz) 

2 2 + zz’ 


(3.30) 
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and from the (.. ,) 3 z> 3 2 component we find 


H’(zz) 


±N ( p - 4 H(zz) 

6 2 + zz ’ 


which gives the nontrivial constraint 


(3.31) 


= =f2, 


(3.32) 


verifying our prediction from more general EFTs (3.25), as well as a functional form for H , 


H 


1 

1 + zz /2 


(3.33) 


In terms of spherical coordinates and standard generators, the ansatz is (again only valid in 
the upper patch): 


A = r 3 (c/)(r)dt ± (1 — cos 9)d(p) 

+ w(r) [(cos (f)d9 — sin 9 sin fidq i>) r 1 =F (sin (j)d6 + sin 9 cos (ftdcft) r 2 ] (3.34) 


Viewed from above, this corresponds to an order parameter configuration which is pointing 
in the same direction across the entire patch. 

With this ansatz one can work out the background equations of motion as a set of ODEs. 
They are similar to those in planar coordinates with additional terms coming from the radial 
slicing curvature and the magnetic field. Solutions with spontaneous symmetry breaking still 
exist below a critical temperature and at large enough chemical potential. Note that due to 
the magnetic field at low temperatures on the sphere the system will develop a AdS 2 x S 2 
near horizon geometry instead of a Lifshitz throat. 


4 Transport 

To study transport properties we look at perturbations on top of the background discussed 
in the previous section. We will only consider hnite frequency but zero spatial momentum 
perturbations of the form 54>(r) e~ lujt . This will allow us to separate the scalar, vector and 
tensor sectors following the mode decomposition of [17]. For the sake of brevity, we restrict 
to the case of p + ip. Since time reversal switches the chirality of the VEV, the results for 
p — ip can be readily derived from what follows by taking uj —» —c u and taking the complex 
conjugate of the fields. 

4.1 Viscosity and tensor perturbations 

In this section we calculate the shear and Hall viscosity of our system (conformal symmetry 
causes the bulk viscosity to vanish). To caclulate the stress viscosity transport coefficients, 


10 




we look at metric perturbations which give us boundary stress response [18]. More explicitly, 
metric fluctuations behave asymptotically as 




—iujt 


+ • • • + h yu e 


—iut 


/r + ... 


(4.1) 


Here we have picked the gauge Sg rr = 5g rfl = 0, and so g runs over t, x, y. When 7 p(r is the 
only nonnormalizable component that is nonzero we have the stress response, 

(T^)W> = (4-2) 

To measure the Hall and shear viscosities, consider a tensor perturbation where the only 
nonnormalizable perturbation turned on is 7 xy . 

(- T xx {-u)T xy (u)) = - iujrm = ( 4 - 3 ) 

(T xy (-uj)T xy (uj)) = - iurj sh = ~ 2 ( 4 - 4 ) 

I xy 

How do such perturbations behave in the superfluid phase? The tensor modes on the p x + ip y 
background are 


$gxy ■ Or, 5A y -\-5A x . o 2 1 &gxx &9yy ■ ® 3 j &A X 5A y . 04 . 

The linearized equations take the form 

d 2 a + K ■ d r a + M ■ a = 0, 


where 
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( Ad ^ 0 0 \ 

9 Z 

k 2 0 0 

0 0 K\ 

0 0 -j£ k 2 
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0 

4 za 2 cj ' u ; 

g 2 fh 

g 2 f h 

m 2 

M 3 

iujwcf) 
r 2 fh 

2iuj(p 

fh 
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4 ihL 2 uw(j> 

M x 

4 K Z W(f)* 

g 2 fh 

g 2 fh 

iuw(j) 
\ r 2 fh 

2iu)(p 

fh 

m 2 
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(4.5) 

(4.6) 


(4.7) 
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fh r 2 f fh' 


Consider the solution which behaves asymptotically as: 

9 1 A 9 Ol 3 0 2 1 

0\ = OioV + —A CJ Oio H-h ..., o 2 = -h 

2 r r 

The final formulae for the viscosities are 


033 


CZ 41 


, a 3 — ——(-..., 04 — ——h 


Vh(w) 

Vsh(u) 


—3i a 33 

2ook 2 L 2 aio ’ 
—3 i cl \3 

2cuk 2 L 2 aio 


(4.8) 


..(4.9) 

(4.10) 

(4.11) 
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Figure 5: (a) Plot of . (b) Plot of u ^. 


4.1.1 Hall viscosity 

We first look at the Hall viscosity, which we compute by numerically integrating the equations 
of motion and using the formulae (4.11). The Hall viscosity vs temperature results are 
given in plot 5a. We see that its magnitude increases as we approach zero temperature. 
Normalizing Hall viscosity value by the superfluid density (figure 5b), we see that the ratio 
equal to —1/2 (±10 - ', within the numerical error we work at) over the entire temperature 
range and for all values of the coupling constant. This is the behavior that was predicted at 
zero tempereature from an effective field theory (2.5). Here we see that its validity is more 
general and implies that at finite temperature, there is no contribution to Hall viscosity 
coming from the normal fluid sector. This is not surprising, as our normal fluid does not 
break parity (see [19] for a holographic model of a normal relativistic fluid which exhibits 
Hall viscosity). We do not have an analytic proof that r] H /p s = —1/2 but suspect that it 
should be possible, likely in a manner similar to the proof of the shear viscosity bound [20]. 

A possible way of understanding this ratio is to note that the ratio of the Hall viscosity to 
the “orbital angular momentum per particle” has been shown to be equal to 1/2 in gapped 
or topological phases at zero temperature [21]. However, if we consider the fact that the 
vector-like order parameter carries angular momentum equal to one, then the superfluid 
density would be equal to the orbital angular momentum per particle and our numerical 
result can be considered as a generalization of the viscosity-angular momentum relationship. 

4.1.2 Shear viscosity 

The plots for shear viscosity vs. temperature at various values of the coupling constant can 
be seen in figure 6a. As expected this dissipative effect vanishes in the zero temperature 
limit. Near T c , where the background geometry is Reissner-Nordstrom, we know that the 
ratio r] s h/S = 1/47T. The plot of this ratio at other temperatures is given in figure 6b. 
We see that as the temperature is reduced, this ratio increases, satisfying the bound [20]. 
However, very close to the zero temperature limit the ratio starts to decrease again. We do 
not understand the origin of this dip. However, we can deduce that it must be coming from 
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Figure 6: (a) Plot of shear viscosity normalized by the chemical potential, (b) Plot of 47Tr> s SH . 
Note that the dip near zero temperature is also present at other values of g. 

subleading terms as the shear viscosity and entropy must have the same scaling behavior in 
the IR Lifshitz fixed point. 

4.2 Conductivity and vector perturbations 

We consider the six following vector modes with harmonic time dependence: 

Sgtz, Sg ty , SA Z X , SA 3 y , SA), SA 2 t . (4.12) 

where we have gauged away the other vector modes. We are interested in calculating the 
electric conductivity, defined by 

Ji = VijE 1 . (4.13) 

It is useful to work in the chiral/antichiral basis as these modes decouple from eachother. 

x + iy + r 1 ± ir 2 

21 = -7^-j T = -7=-, 

V2 y/2 

This allows us to directly compute the eigenvalues of the conductivity matrix, 


®xx i xyt ®xx I'&xy 


(4.14) 

Our bulk modes are similarly organized as 



dA *^^, SA± = dAlt ^ A \ Sg tz = 

$Qtx i'&Qty 

V2 ' 

(4.15) 

Dehning 

a= ( SA 3 Z , «54 t + , Sg tz ) a= ( SA |, SA~, 

Sgtz ), 

(4.16) 
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the a and a modes are completely decoupled. In what follows, we restrict the discussion to 
the a modes. The linearized equations for the a modes are given by 4 : 

d 2 r a x + K-d r a + M-a = 0, (4.17) 


where 


where 
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g 2 r 2 ((f)—Uj) 
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uj 2 2 K 2 ww'(j)' 

fh g* r 2 ( 0 _ w) _^ 2 / 


2 k 2 0' 2 / 2 K 2 cj)' 2 (j) + 

g 2 h ’ m V^ 2 (0-(n) + / J ' 


At large radius, the fields (4.16) will be of the form 


_ , / , _ , a 30«h 

0,1 — a io + On r + • ■ ■, 02 — 020 H- 
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anWoUioWi 
L 2 (/i — a;)r 


(a 2 o — a3QWo)p 
(/i — cj)r 


2 , ( , wo(anWo - aio^i) , wop(a, 2 o 

03 — a 30 r + I a 10 p H- — -7- 1 - 

V w(u-co) u - 


■ n.onllln'W 2 K 2 


(4.18) 


It is important with these perturbations do not only turn on a linearized non-normalizable 
Al 2 by way of aio but also possibly Af and g tz . These may be removed by a linearized gauge 
transformation [11] and diffeomorphism and do not change the equations of motion. However, 
they affect the formulae for conductivity. Consider a linearized gauge transformation followed 
by a diffeomorphism 

A = e _<wt A + t “, £ = e~ iut i z d- z , (4.19) 


which we fix by requiring after the gauge transformation A+, g tz r^i o(i/ 0 - 


This is done 


via 

tz * a 30 , + i (°2o~~ 2a 30 u’o) 

S • ^ Cl 

co 2 ii — oo 

After this fixing we find the conductivity is 5 


(4.20) 


a z _ 1 On{p — u>) + (<220 — ^o^qWo)wi (4 21 ) 

a n iu ai 0 (/i — a;) + w 0 (a 2 o — 2a 30 w 0 ) ’ 

4 The a mode equations are simply derived by a time reversal transformation which takes uj —> — w and 
conjugates the fields (i.e. it takes a —> a). 

5 It is clear from this formula that there will be a pole in the imaginary part of the conductivity at uj = g. 
This is isospin precession which is undamped and temperature independent and follows directly from the 
nonabelian Ward identities [22, 23]. 
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(4.22) 


where we have normalized by the value of the conductivity at co —> oo, 

, ,_L 2 

@n ®xx y t OO) . 

r 

Since we are interested in backgrounds with spontaneous symmetry breaking where Wq = 
0, these formulae simplify and we have: 




an 


+ - 


U7ia 2 o 


(4.23) 


<r n icoaio iojci io(h — ta) 

The conductivity for the anti-chiral modes is derived from this by complex conjugation and 
taking u> —> —u: 


Oj_ 

a, 


a ii W1O-20 


n 


icoa w icoa w (/d + uj) 


(4.24) 




Figure 7: Plots of yy (left) and yy (right) for various values of T/T C - 


The typical plots for the spectral densities can be seen in figure 7. We note that these 
are qualitatively the same as their strict probe counterparts 6 in [11], We see that there 
is a pseudogap, i.e. a depletion of states as we approach zero temperature, which is a 
characteristic of holographic superconductors. The respective plots for the conductivities 
can be seen in figure 8. Near T c , where the order parameter vanishes, parity is restored, 
which implies yy and \z are equal. This implies that the DC conductivity must vanish. 

By comparing the values of o xy to the superfluid density in the zero temperature limit, 
we can read off the coefficient of the non-universal £ term in the effective held theory (eq. 
(2.7)). In the probe limit, this ratio is exactly equal to one. Turning on the back reaction of 
the geometry, we see that for large values of g, the deviation from probe goes as g~ 2 . A fit 
to the ratio for various g reveals that Co = 1 — 1.15 g~ 2 + 1.05 g~ i + 0(g~ 6 ). As expected, 
this coefficient is non-universal and depends on the details of the theory and can also vary 
by temperature. We suspect that this is the reason why unlike the ratio of Hall viscosity 
to superfluid density which is independent of temperature, the same cannot be said for Hall 
conductivity. 

6 The one exception is the dip in the spectral density near zero frequency at Tq- In the probe limit the 
dip vanishes. 
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Figure 8: Plots of conductivities at zero temperature. The orange solid line is the real part 
and the blue dashed line denotes the imaginary part. 


5 Conclusions and outlook 

In this paper we analyzed a chiral superfluid model and computed response functions away 
from the probe limit. While some characteristics of the system (e.g. shape of electric 
conductivity plots) remain the same, there are some quantities that can only be calculated 
away from the probe limit. Due to the emergent Lifshitz geometry that governs the infrared 
dynamics, many of the thermodynamic variables are quantitatively different. 

Looking at the viscosity response of the system, we found that the Hall viscosity to 
superfluid density ratio matches the zero temperature EFT prediction even at finite tem¬ 
perature, as does the superfluid shift. We also computed the shear viscosity and found that 
the viscosity-entropy ratio respects the holographic bound, but has a feature near zero tem¬ 
perature where it decreases. Our results for transport are numerical only and it would be 
interesting to understand them from a first principle approach. 

Unlike the low energy EFT analysis which only kept track of the Goldstone degree of 
freedom, this holographic model has the capacity to describe vortices. While these have been 
studied in the probe limit near critical magnetic fields [24] it would be very interesting to 
study them with back reaction. It would also be interesting to study the structure of bound 
states that arise on different topologies either due to the vortices created from a mismatch 
in the magnetic flux or in the presence of boundaries. 
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